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Abstract 

Let p the density matrix of a mixed Gaussian state. Assuming 
that one of the Robertson-Schrodinger uncertainty inequalities is sat- 
urated by p, e.g. (A' Xxf (A? Pi) 2 = A^X^Pi) 2 + ±/i 2 , we show that 
there exists a unique pure Gaussian state whose Wigner distribution is 
dominated by that of p and having the same variances and covariance 
A p Xi, A p Pi, and A p (Xi,Pi) as p. This property can be viewed as 
an analytic version of Gromov's non-squeezing theorem in the linear 
case, which implies that the intersection of a symplectic ball by a single 
plane of conjugate coordinates determines the radius of this ball. 



1 Statement of Results 

Consider a mixed quantum state, identified with its density matrix p, and 
let 

Wp(x,p) = (^) n / e-f^ix + \y\p\x - \y)dy 



be its Wigner function. Viewing the latter as a quasi-probability distribution 
we assume that it satisfies 

(1 + \z\ 2 )\W p{z)\dz < oo (1) 
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(we have set z = (x,p)) so that the first and second moments of Wp exist. 
The covariance matrix (CM) of p is defined by 

E p = [ (z-(z))(z-{z)) T Wp(z)dz (2) 

where (z) = Tr(zp) is the mean value vector. It will be convenient to write 

_ (AP(X,X) AP(X,P)\ 
p \AP(X,P) AP(P,P)J {) 

where AP(X,X) = (A"(X i; Xj))^^, AP(X,P) = (AP(X h ^))i<vj<»> 
A P (P,P) = (A p (Pi, Pj))i<i,j<n where A p (X.i, Xj), etc., are the covariances. 
Setting (APXj) 2 = A p (Xj',Xj) and (APPj) 2 = AP(Pj,Pj) the Robertson- 
Schrodinger (RS) inequalities 

{APXjfiAPP,) 2 > APiXj,^) 2 + \h 2 (4) 

hold for j = 1, n. We will say that these inequalities are partially saturated 
if at least one (but not all) of them are equalities. 

We now assume that the state p is Gaussian; this means that 



1 v n/2 

Wp(z)={ — \ (detsr^exp 



\{z-(z)) T Tr\z-(z)) 



(5) 



(thus £ is assumed to be invertible, but this is no restriction: see section [2]). 
We will show that: 

Theorem 1 Suppose that anyone of the (RS) inequalities is saturated, 
for instance 

(APX 1 ) 2 (APP 1 ) 2 = AP(X 1 ,P 1 ) 2 + lh 2 . (6) 

Then: 

(i) There exists a unique pure Gaussian state such that 
A^X X = APX X , A^Pi = A P P X , A^(Xi,Pi) = A p (Xx,Pi) 

(ii) That state ip is the only Gaussian whose Wigner function satisfies 
the inequality. 

Wip(x,p) <Wp(x,p). (7) 
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We emphasize that no assumptions are made on the other variances or 
covariances of the state p: they may for instance unknown. 

One of the crucial steps in the proof is of a topological nature, and is 
deeply related to Gromov's [16J symplectic non-squeezing theorem (nick- 
named by V.I. Arnol'd [3] the "principle of the symplectic camel" ; see [9j 
for a discussion of that terminology). Let us now call "symplectic ball" the 
image S(Br) of a phase space ball Br with radius R by a linear (or affine) 
symplectic transformation S; the number R is the radius of the symplectic 
ball. Since symplectic transformations are volume preserving, the volume 
of S(Br) is the same as that of the ball Br. More surprising is the fol- 
lowing rather unexpected property: the section of the ellipsoid S(Br) by 
any plane of conjugate coordinates Xj,pj passing through its center is al- 
ways irR 2 . Thus, a single "tomography" of a symplectic ball unambiguously 
determines its radius! This is counterintuitive, because one would expect 
that the area of the section of an ellipsoid by different planes yields different 
results. It turns out that this property (of which we give a proof in the 
Appendix) is equivalent to the linear version of Gromov's theorem (see [15] 
for a review of Gromov's theory; also [23] for various developments). Its 
relation with Theorem [1] comes from the fact that the Wigner transform of 
a Gaussian is a Gaussian of a very special type: its covariance ellipsoid is 
a symplectic ball, and the state is thus entirely determined by a section of 
this ellipsoid by a single plane of conjugate variables. 

This paper has two forerunners [5] and [6] where we showed that it was 
possible to associate a unique pure Gaussian state to a mixed state satisfying 
a certain topological condition related to the uncertainty principle, see (fl~2j) 
in section [2j 

Notation. We will use the following notation and terminology: The generic 
variable of phase space E? n is z = (x,p) with x = (x±, ...,x n ), p = (pi, ...,p n ); 
in calculations x,p,z will always be viewed as column vectors. Sp(2n,M) is 
the symplectic group; it consists of all real 2n x 2n matrices S such that 



Dirac's notation h for h/2ir (h Planck's constant). 

Acknowledgements. The present work has been supported by the Aus- 
trian Research Agency FWF (Projektnummer P20442-N13). 



S T JS = J (or, equivalently, SJS 



•r 



J) where 




We use 
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2 Geometry and Uncertainty 



In his seminal paper "Geometry and Uncertainty" [21] F. Narcowich (also 
see [22\ 126] ) proves, among other things, the following deep result: a real 
symmetric matrix S is the covariance matrix of a quantum state if and only 
if 

ih 

£ + —J is positive semi- definite (8) 

(abbreviated as £ + |j > 0). Condition © implies ([2~T]. Lemma 2.3) that 
£ is fact definite positive, and hence invertible. It will be convenient to use 
the auxiliary matrix 

M = -S- 1 (9) 



in which case (ED becomes 



M _1 +iJ>0. (10) 



Narcowich's result seems to have become "folk wisdom" (and is therefore 
not very much cited) ; its proof uses arguments from hard harmonic analysis 
(a symplectic version of Bochner's theorem on positive measures [171 119|, [20] ) 
and is often stated (most of the time without proof) in the quantum-optical 
literature; see for instance the references in pQ. A caveat: as we have shown 
in |13) there exist self-adjoint operators with trace one whose covariance 
matrix satisfies condition ([8]), but which are not non-negative, and hence 
do not represent a quantum state. However, given a covariance matrix E 
satisfying (|8|) one can always construct a quantum state with covariance 
matrix £, namely the Gaussian state p whose Wigner transform is 

Wp(z) = (det M) exp (^-^z T Mz 

In [3 El El EE HO] we have shown that condition (E|) is equivalent to a 
topological statement involving the symplectic capacity of the covariance 
ellipsoid 

n = {z: \z T YT x z < 1} = {z : z T Mz < H}. (11) 

In fact, the algebraic condition M^ 1 +iJ > 0. is equivalent to the following 
property: 

The symplectic capacity c(f2) of the covariance ellipsoid is > ^h. (12) 

The symplectic capacity c(f2) is defined as follows: if n = 1 it is just the 
area of the ellipse f2; in higher dimensions it is the supremum of all numbers 



4 



7rR 2 such that a symplectic ball S(Br) is contained in Q. In [12] we have 
called a symplectic ball with radius \f% a "quantum blob" . 

Both conditions (|5J), (fl2j) are, for a given quantum system, equivalent to 
the Robertson-Schrodinger (RS) inequalities 

(WXjf^Pj) 2 > A^Xj,^) 2 + \h 2 , 1 < j < n. (13) 

Using either © or (JT2J) it is easy to show that the RS inequalities are co- 
variant {i.e. they retain they form) under linear or affine symplectic trans- 
formations. If we set (X',P') = S(X,P) where S is a symplectic matrix, 
then 

(A^) 2 (A^Pj) 2 > A"(Xj,Pj) 2 + \h 2 , 1 < j < n (14) 
since the new covariance matrix is £' = ST I S T and 

because SJS T J = J as 5 is symplectic. 

The key to the proof of the equivalences (|8|)^=> (IT2]) ^=> (jT3"]) is the fol- 
lowing well-known symplectic diagonalization result (Williamson's |25| the- 
orem, of which |24| gives an elementary proof). Let M be a (real) symmetric 
positive-definite matrix of size In. There exists S G Sp(2n,R) and positive 
numbers Ai, A n such that 

S T MS=(^ ^ , A = diag(Ai,...,A n ). (15) 

The numbers Ai, A n are called the symplectic eigenvalues of M (or some- 
times Williamson invariants); they are written in decreasing order: Ai > 
• • • > A n and the array 

Spe C(T (M) = (Ai,...,A n ) (16) 

is the symplectic spectrum of M. The symplectic eigenvalues are calculated 
as follows: consider the product JM of M by the standard symplectic matrix 
and let M 1 / 2 be the positive square root of M\ the matrices M 1 / 2 JM 1 ! 2 and 
JM are equivalent, and hence have the same eigenvalues. Since 
is antisymmetric (because J is), these eigenvalues occur in pairs (iXj, ±^Aj), 
Xj > 0, and the \j are precisely the symplectic eigenvalues of M. 

The symplectic spectrum has the following straightforward properties: 
if c> 

Spe C(J (cM) = (cAi,...,cA n ) (17) 
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and 

Spec ff (M- 1 ) = (A; 1 ,...,Ar 1 ). (18) 

The relation between symplectic spectrum and symplectic capacity is essen- 
tial: let M be a symmetric positive definite matrix and consider the phase 
space ellipsoid Q : z T Mz < H. We have 

c(0) = ^ if Spec CT (M) = (Ai,...,A n ). (19) 
Al 

3 The Wigner function of a Gaussian 

We will consider (normalized) Gaussians of the type 

iI) X ,y{x) = (TTh)- n/i (detX) 1 ^e-^ xTix+iY)x (20) 

where X and Y are real symmetric n x n matrices, X positive definite. 
We will need the two following properties of the Wigner function 

W^(x,p) = (^) n / e-i^(x + \yW{x - \y)dy (21) 

(for a complete review of the Wigner formalism we refer to Littlejohn [18] 
and de Gosson [11] ) : 

Symplectic covariance. Let S £ Sp(2re,M) and S be any of the two 
metaplectic operators associated with S. Then 

W(Sil>)(z) = W^{S- l z). (22) 

Recall that the metaplectic operators are defined as follows: the symplectic 
group Sp(2n,R) has a covering group of order two, the metaplectic group 
Mp(2n,R). That group consists of unitary operators on L 2 (IR ra ), and is 
generated by the following elementary unitary transformations: 

• The modified Fourier transform 

(with the convention i 1 / 2 = e l7r / 4 ); 

• " Chirps 1 ' 

Vpip{x) = exp ( -^-x T Px J , P = P T ; 
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• Dilations 

M L %l){x) = Vdet Lip(Lx) , detL/O 

where V det L = i m ^/\ det L\, m = or 2 if det L > and m = 1 or 3 
if det L < 0. 

For a detailed study of the metaplectic group Mp(2n,M) see [7]. 
Transformation of Gaussians. The Wigner transform of a Gaussian is 
itself a Gaussian. In fact 

WiI> x ,y( z ) = (vr^)-"e-^ TGz (23) 
where G is the real 2n x 2n matrix 

It turns out that the matrix G defined in the formula above is both 
positive-definite and symplectic: we have G = S T S where 

S= (x-^Y X' 1 / 2 ) (25) 

and S is obviously in Sp(2n,M). The argument can be reversed: given a 
positive-definite symplectic matrix G, we can always find S G Sp(2n,IR) 
such that G = S T S where S has the form (|25|) : this determines X and Y 
and hence a Gaussian state tpx,Y satisfying (f23|) . 



4 Proof of Theorem [T] 

It is sufficient to assume that p is centered at the origin, that is (z) = 0. 
First step. We begin by noting that the saturation of one of the RS in- 
equalities, for instance 

(APX^iAep,) 2 = dJ>{X u Fx) 2 + \h\ (26) 

implies that c(O) = \h. Let us prove this by reductio ad absurdum. Suppose 
that c(S7) > and let D = S T MS be a symplectic diagonalization (fT5|) 
of M = Since c(O) is a symplectic invariant {i.e. c(5(0)) = c(f2) 

for every S G Sp(2n,M)), our assumption can be rewritten c(A) > 
where A is the ellipsoid defined by z T Dz < h. In view of formula f)19f) 
we have c(A) = irh/Xi, recalling that the Xj form a decreasing sequence. 
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The assumption c(A) > \h is thus equivalent to Ai < 1 and hence to \j < 1 
for all j = 1, n. Consider now the matrix 



B~ x + iJ 



A- 1 I 
-I A- 1 



> 0. 



Its eigenvalues are the roots of the polynomial 



and are thus the real numbers tj = A 1 ± 1. Since Xj < 1 for every j this 
means that tj > for every j and hence D^ 1 + i J > 0; returning to the 
covariance matrix, it follows that we also have £ + y J > 0. But this leads 
to a contradiction: in view of Sylvester's criterion for definite positiveness 
the condition £ + y J > implies the following condition on the minors of 
S + f J: 



ill 
2 



> 



that is 



(A^) 2 (A^) 2 > A^A^Pj) 2 + \h 2 



for all j = 1, n. 

Second step. Let us show that if Q is a phase space ellipsoid such that 
c(f2) = |/t, then $7 contains a unique quantum blob (i.e. a symplectic ball 
with radius Vh); we do not assume explicitly that Q is a covariance ellipsoid: 
the argument is quite general. We begin by noting that the condition c(fi) = 
\h implies, by definition of a symplectic capacity, that $7 contains the image 
by some S G Sp(2n, M.) of a ball (and no symplectic ball with larger 
radius). This symplectic ball can be explicitly constructed: assume that f2 
is given by z T Mz < h, M > 0. In view of Williamson's diagonalization 
theorem, there exists S G Sp(2n, R) such that 



S T MS = D 



A = Spe C(T (M) = (Ai,...,A n ; 



(27) 



(cf. (|15p . The inequality z T Mz < h is equivalent to 



^2\j(xf +pf) < h , (x',p) = S{x,p) 



(28) 
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and we have c(O) = nh/Xi = -^h so that Ai = 1; it follows that the el- 
lipsoid (f28|) contains B^ hence Q contains the symplectic ball S(B^). 
There remains to prove the uniqueness of the a quantum blob contained in 
Q. Using if necessary a phase space translation it is sufficient to consider 
the case where B/^ : \z\ < H (in view of elementary geometric considera- 
tions, the largest symplectic ball contained in Q must be centered at the 
origin: see [7j, §8.4). Let us assume that there exists S' £ Sp(2re,M) such 
that S'(B^f t ) is also contained in the ellipsoid ft; we are going to show that 
S'(B^) = S(B^qi), following the argument in de Gosson P2EHI5J. Conju- 
gating if necessary S and S' with an adequately chosen symplectic matrix 
we may assume, using Williamson's theorem, that the Xj are the symplectic 
eigenvalues of the matrix of ft. The condition c(f2) = then means that 
Ai = 1 in view of formula (fl~9|) . Let us show that U = S'S -1 belongs to the 
group of symplectic rotations 

U(n) = Sp(2n,M)nO(2n,M); (29) 

the claim will follow since then U(B^) = B ^ so that S(B^) = S'(B^). 
Clearly U £ Sp(2n,M) so it suffices to show that U is in addition a rotation; 
for this it is sufficient to check that UJ = JU. Set R = D 1 I 2 UD~ 1 / 2 
where D = diag(Ai, A n ); we have U T DU = D hence R T R = I so that 
R is orthogonal. Let us prove that R is in addition symplectic. Since J 
commutes with every power of the diagonal matrix D we have, taking into 
account the relation JU = (U T )~ 1 J (because U is symplectic): 

JR = D^JUD- 1 ' 2 = D l l 2 (U T y x JD- l l' 2 
= D 1 l 2 (U T )- l D- x ' 2 J = {R T )~ l J 

hence R T JR = J so that R is in Sp(2n,R). Since R is also a rotation we 
have R 6 U(n) and thus JR = RJ. Since U = D~ l / 2 RD 1 / 2 we have 

JU = JD~ 1/2 RD 1/2 = D~ 1/2 JRD 1/2 
= D~ 1/2 RJD 1/2 = D~ 1/2 RD 1/2 J 
= UJ 



proving our claim. 

Third step. We now define a normalized Gaussian state ip by specifying 
its Wigner transform: 



z T {SS 1 )- L z 



T\-l, 
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where S € Sp(2n,M) is defined as above, i.e. S{B^} is the largest quantum 
blob contained in Q = {z : z T Mz < h}. We note that the choice of S 
is irrelevant, because we have seen above that if S(B^) = S'(B^) then 
S' = SU with U e U{n) so that S'S' T = SS T . Since S(B^ C we have 
z T Mz < z T (SS T )~ 1 z for all z G R 2n and hence, taking definition (J5j) of 
Wp(z) into account, 

Wtp(z) < Wp(z) (30) 

which is the inequality ([7]) in Theorem [TJ That ip is the only Gaussian 
state satisfying this inequality follows at once from the uniqueness of the 
quantum blob S(B^). Let us finally show that the states ip and p have 
the same covariances (A p Xi) 2 , (A p Pi) 2 , and A p (Xi, Pi). The covariance 
matrix of ip is £^ = ^SS T and the inclusion S(B^) C ft implies that 
we have 0* C S(B^)* where $7* and S(B^)* are the dual ellipsoids of 
and S(B^), obtained by using a Legendre transformation (for a review of 
the latter see [2Y]). Since the dual ellipsoid of an ellipsoid ^z T Y,z < 1 is 
\z t Yj~ x z < 1 this implies that S(B^)* = (S T )~ l B 2 ,^, so the points of 
these dual ellipsoids satisfy 

j Z T SS T z < ^z T ^ p z < 1. (31) 

Let us cut = {S T Y 1 B 2/VTi and Q* = {z : ^z T J: p z < 1} by the 

plane P^i of coordinates x\,p\. This section consists of two ellipses: P\ t \ n 
{S T )~ l B 2 i^j i , which has area An/h, and 

(A^) 2 A'p^P^N /xA 

a^x^px) (A^po 2 ; \ Pl J -> 

taking the saturation of the first RS inequality ([6]) into account, P^i n Q* 
also has area 

2vr [(A' , X 1 ) 2 (A^ > P 1 ) 2 - A'C^Pi) 2 ] -12 = Att/H. 

Since P^i n ST C P^i n S(B^)* equality of areas implies that the ellipses 
Pi t i n ri* and Pi t i n S{B r^)* are identical, hence the covariances (A p Xl) 2 , 
(A p Pi) 2 , and A p (Xl,Pi) are the same for both states ip and p. QED. 

Appendix: The Linear Gromov Theorem 

Let us prove the following linear version of Gromov's non-squeezing theorem. 
We are following de Gosson [IT], §5.1.2; we give two independent proofs. 
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Proposition 2 Let S £ Sp(2n,M) and Bp : \z\ < R The intersection of 
S(Bfi) by a plane of conjugate coordinates Xj,Pj is an ellipse with area irR 2 . 

First proof. It relies on the fact that the form pdx = ^ZjPjdxj is a relative 
symplectic integral invariant, that is: if <p is a symplectomorphism of M? n 
and 7 a loop in M 2n , then 



pdx = d> pdx (32) 

7 J <K7) 

(see for instance Arnol'd [2j, §44, p. 239). We claim that the ellipse Yj = 
S(Bp) n Vj, intersection of the ellipsoid S(Bp) with any plane Vj of con- 
jugate coordinates Xj,pj has area ttR 2 ; the proposition immediately follows 
from this property. Let jj be the curve bounding the ellipse Tj and orient 
it positively; the area it encloses is 

Area(rj) = <* pdx = f pdx = ttR 2 (33) 

H Js-Hfi) 

(because S ,_1 (7 J ) is a big circle of Bp); notice that the assumption that Vj is 
a plane of conjugate coordinates essential for (j33|) to hold, making 

the use of the formula (|32|) possible [more generally, the argument works 
when Vj is replaced by any symplectic plane] . 

Second proof. With the same notation as above we note that the set 

s- 1 [S(Bp)nv j ] 

is a big circle of Bp, and hence encloses a surface with area irR 2 . Now, Vj 
is a symplectic space when equipped with the two-form Uj = dpj A dxj and 
the restriction of S to Vj is a linear symplectomorphism from (Vj,o~j) to the 
symplectic plane S(Vj) equipped with the restriction of the symplectic form 
a. Symplectomorphisms being volume (here: area) preserving, it follows 
that S(Bp) n Vj also has area ttR 2 . 

Remark. It would certainly be interesting to generalize the first proof to 
arbitrary symplectomorphisms. The difficulty comes from the following fact: 
the key to the proof in the linear case is the fact that we were able to derive 
the equality 

/ Pjdxj = ttR 2 

by exploiting the fact the inverse image of the Xj,pj plane by S was a 
plane cutting Bp along a big circle, which thus encloses an area equal to 
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irR 2 . When one replaces the linear transformation S by a non-linear one, 
the inverse image of the Xj,pj plane will not generally be a plane, but 
rather a 2-dimensional symplectic manifold. If the following property holds: 
The section of Br by any 2-dimensional symplectic manifold containing the 
center of Br has an area at least irR 2 then we would have, by the same 
argument as above 

Pjdxj > ttR 2 . 

We do not know any proof of this property; nor do we know whether it is 
true! 
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